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Operation of nanomechanical devices in water environment has been challenging due to the strong 
viscous damping that greatly impedes the mechanical motion. Here we demonstrate an optome¬ 
chanical micro-wheel resonator integrated in microfluidic system that supports low-loss optical reso¬ 
nances at near-visible wavelength with quality factor up to 1.5 million. The device can be operated 
in self-oscillation mode in air with low threshold power of 45 uW. The very high optical Q allows 
the observation of the thermal Brownian motion of the mechanical mode in both air and water 
environment with high signal-to-background ratio. A numerical model is developed to calculate the 
hydrodynamic effect on the device due to the surrounding water, which agrees well with the exper¬ 
imental results. With its very high resonance frequency (170 MHz) and small loaded mass (75 pg), 
the present device is estimated to have mass sensitivity of attogram level in liquid environment with 
bandwidth of 1 Hz. 


I. INTRODUCTION 

Micro- and nano-scale mechanical resonators have been 
developed as important tools in both fundamental stud¬ 
ies and technological applications. Because of their very 
small spring constant and mass, even tiny forces acting 
on the resonator or masses adhering on its surface can 
greatly alter its dynamics, which can be detected with 
extremely high precision. Force sensitivity down to sub- 
attonewton (< 10 -18 N) level has been demonstrated [T], 
which allows detection of the force from a single electron 
spin [5] . As the dimensions of the mechanical devices con¬ 
tinue to be scaled down, the mass sensitivity has reached 
attogram (10 -21 kg) p] 0], zeptogram (1CT 24 kg) [5], 
yoctogram (10 -27 kg) [B| level, and down to the mass 
of a single proton [7j in the past few years. A recent 
demonstration showed that this nanomechanical technol¬ 
ogy holds a promise to perform mass spectrometry on a 
single molecule with extremely high resolution [5]. 

However, so far most of the high performance nanome¬ 
chanical sensing systems operate in vacuum or air, only 
very few operate in liquid due to the rapid deteriora¬ 
tion of device performance caused by the enormous flu¬ 
idic damping. This dramatically lowers the quality factor 
(' Q ) and the fluid displaced by the movement adds mass. 
Both effects make operation in liquids extremely chal¬ 
lenging. Moreover, the fluidic damping becomes increas¬ 
ingly significant when the resonator dimension is scaled 
to micro- and nano-scales since the viscous force scales 
inversely to the square of the characteristic length, i.e., 
/rV 2 u ~ L~ 2 . For example, in tapping mode AFM the 
typical quality factor in liquid are low: Q < 5 mm- 

To circumvent the severe damping problem, one ap¬ 
proach is to use a hollow resonator structure where fluids 
and target analytes flow in a fluidic channel embedded 
inside mm- The resonator can then be operated in 
vacuum or ambient air so that the mechanical damp¬ 
ing is not compromised. With this approach attogram 
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sensing has been achieved [13]. Nevertheless, operation 
of nanomechanical system in a genuine liquid environ¬ 
ment remains important since only then the sensing can 
take place in environments where many biological and 
chemical samples naturally reside in M- Efforts towards 
this direction include the development of transduction 
schemes that can efficiently actuate and readout the mo¬ 
tion of the resonator in the highly dissipative liquid envi¬ 
ronment. Transduction methods such as thermo-optical 
excitation mm, magnetomotive drive and detection 
m, and piezoelectric actuation |TB] ) have been demon¬ 
strated. On the other hand, a recent demonstration sug¬ 
gests that the use of higher modes of a cantilever could 
be a promising candidate for operation in liquid as they 
generally have higher quality factors and a smaller mass¬ 
loading effect from the water [19] . 

Besides sensing applications, efficient transduction of 
nanomechanical resonators in viscous environment also 
benefits the understanding of fluid dynamics at small 
length scales. Nanomechanical resonators have been an 
important tool for exploring new parameter regimes in 
fluid dynamics, such as the high frequency regime |20l 
m, transitions between Newtonian and non-Newtonian 
regimes [20l 122] and between hydrodynamic and kinetic 
regimes [23]. It also enables the study of stochastic dy¬ 
namics of fluid-structure interaction due to Brownian 
noise which can be used for developing highly 

sensitive AFM in liquid 1211- 

In this paper, we develop an optomechanical micro¬ 
wheel resonator integrated in microfluidic system that 
operates efficiently in fluidic environment. The resonator 
supports high Q optical whispering gallery inodes, which 
transduce the mechanical motion with ultra-high sen¬ 
sitivity and allow optomechanical self-oscillation in air 
with low threshold power. Using the highly sensitive 
cavity-enhanced optical readout, we are able to resolve 
the thermomechanical noise of the mechanical radial 
breathing mode in water with high signal-to-background 
ratio. To explain the changes in resonance frequency and 
damping rate caused by the liquid, a boundary integral 
method is developed to calculate the hydrodynamic func¬ 
tion for axial-symmetric systems. Together with its small 
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effective mass (75 pg with the entrained liquid included), 
high resonance frequency (170 MHz in water), and high 
displacement sensitivity (15 am/VHz), mass sensitivity 
of 2.5 attogram with measurement bandwidth of 1 Hz is 
predicted. This device shows promise in sensing applica¬ 
tions in fluidic environment. Our design of the integrated 
system and the robust fabrication process also provide a 
viable approach for large scale integration of nanopho¬ 
tonic devices with nricrofluidics. 


II. DEVICE DESIGN AND FABRICATION 

One of the technical obstacles of realizing photonic res¬ 
onators in water environment is the optical absorption 
due to the water, which can severely degrade the opti¬ 
cal Q and cause undesirable thermo-optical heating of 
the surrounding water. It is well known that water is 
strongly absorptive in the near-infrared range where most 
integrated optomechanics experiments are done. The op¬ 
tical extinction coefficient at wavelength of A = 1.55 /im 
is k = 9.86 x 10 -5 , [28] which corresponds to a loss of -34 
dB/cm for an electromagnetic plane wave propagating in 
the medium. Therefore, telecommunication wavelength 
commonly used in silicon photonics is not suitable for 
this purpose. Although ID photonic crystal resonators 
at this wavelength range has been demonstrated in wa¬ 
ter, a low quality factor Q ss 3000 was obtained j25]. 
To achieve high Q optical resonances in water, here we 
developed a photonic resonator that operates in near- 
visible wavelength at around 780 urn. The extinction 
coefficient of water at this wavelength is about 700 times 
smaller: k = 1.43 x 10“'. [25] This is equivalent to an 
attenuation of only -0.1 dB/cm, which would result in 
an absorption-limited Q of 10 million. This wavelength 
also falls within the biological transparency window [30] 
and is therefore useful in biology and for medical appli¬ 
cations. For the waveguiding material, stoichiometric sil¬ 
icon nitride is used because of its low optical absorption 
in the visible range, as shown by a recent demonstration 
of optical resonances with Qs up to 3 million in the 652 
- 660 nm wavelength range. m 

The design of the optomechanical resonator in this 
study is a suspended micro-wheel structure [52145H] . Me¬ 
chanically it supports a radial breathing mode, which has 
been demonstrated to have high mechanical quality fac¬ 
tor up to Q m ~ 4000 at 1.3 GHz operated in air [33]. 
The motion of the radial breathing mode does not in¬ 
volve change of center of mass. Therefore, it is expected 
to cause minimum energy loss due to radiation (i.e., mass 
transport) through the fluidic environment. Fig. 0(a) 
shows the displacement profile of the fundamental radial 
breathing mode of a micro-wheel of 10 um radius sim¬ 
ulated using finite element analysis. This mode has an 
effective mass of 65 pg and a resonance frequency of 175 
MHz in vacuum. 

The mechanical mode is optomechanically coupled to 
the optical whispering gallery modes [36]: as the displace¬ 


ment of the mechanical mode vibrates in the radial direc¬ 
tion, it changes the resonance wavelength of the optical 
mode, which can be used as the readout of the mechani¬ 
cal motion. Fig. 0(b) shows a cross-sectional plot of the 
simulated radial electric field distribution of the first four 
transverse-electric (TE) modes in air. The optical mode 
is mainly confined at the outer rim of the micro-wheel. 
Therefore scattering loss due to the presence of the an¬ 
choring spokes at the inner rim has a negligible effect. 
For the case when the device is immersed in water, simu¬ 
lation shows that the fourth TE mode becomes not well 
confined due to the smaller index contrast. 

The device is fabricated in a 200 nm thick stoichiomet¬ 
ric silicon nitride film deposited by low-pressure chemi¬ 
cal vapor deposition (LPCVD). Details of the fabrica¬ 
tion procedures can be found in the supplementary ma¬ 
terial. The device layer is cladded with a 2 um thick 
thermal oxide from the bottom and a 3 um thick plasma- 
enhanced chemical-vapor deposited (PECVD) oxide from 
the top. Microfluidic channels were crafted into the top 
cladding layer directly above the device exposing the res¬ 
onator structure to the external environment. The chip 
was bonded to a thin cover glass slide using PDMS as a 
bonding layer. For each device, an evanescently coupled 
waveguide is fabricated next to the resonator and a pair 
of gratings is used for coupling light into and out of the 
waveguide. Fig. 0 (c) and (d) show SEM images of a 
device with outer radius of 10 um and ring width of 2.5 
um before the cover glass bonding. As can be seen in the 
figure, a tapered structure is designed to provide a ro¬ 
bust support to the free-standing waveguide. In fact, the 
whole structure is sturdy enough that it can be released 
directly in a wet undercut process without the use of 
a critical point dryer. This allows repeatedly switching 
between operation in liquid and gaseous environments. 
Fig. 0 (d) shows a device positioned in the middle of a 
microfluidic channel, which is highlighted in blue color. 
Fig. 0 (e) is an optical image showing the top-view of 
the chip after the glass bonding. A fully packaged de¬ 
vice is shown in Fig. 0 (f), where a fiber array probe 
is approached from the top and aligned to the grating 
couplers on the chip. This design seamlessly integrates 
the microfluidic and photonic systems and allows efficient 
fiber access to the on-chip waveguides. It provides a vi¬ 
able approach for large scale integration of nanophotonic 
devices with microfluidic system. 

The device was characterized using the setup shown in 
Fig. 0 A tunable diode laser (New Focus TLB-6712) 
with a tuning range of 765 - 781 nm was used to probe 
the optical spectral response of the device. A small por¬ 
tion of the laser output was tapped out and sent to a 
wavelength meter for wavelength and intensity calibra¬ 
tion. Using a fiber polarization controller, the laser light 
was adjusted to be TE-polarized to optimize the coupling 
efficiency of the grating couplers. Next, a fiber probe was 
aligned from the top to the grating couplers on the device 
and the transmitted light was collected on a high speed 
photoreceiver (New Focus 1601, 1 GHz bandwidth). Part 
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FIG. 1. (a) Normalized displacement profile of the fundamental radial breathing mode, (b) Electric field (radial component) 

distribution of the first four TE modes (TE;) in air, where the index i represents the mode orders in the radial directions, 
(c) Scanning-electron microscope image showing a top view of the device, (d) Angled view SEM image showing the device 
positioned in the middle of the microfluidic channel (highlighted in blue) (e) Optical top view of the chip showing arrays of 
devices and two microfluidic channels (highlighted in blue), (f) Photo of a fully packaged device aligned to a fiber probe. 



FIG. 2. (a) Schematic of the measurement setup. TDL: 

Tunable diode laser. FPC: Fiber polarization controller. PR: 
Photoreceiver. DAQ: Data acquisition system. 


of the signal was sent to a data acquisition system to mea¬ 
sure the dc transmission while the other part was sent to 
an electrical spectrum analyzer for spectral characteriza¬ 
tion. In this setup, all the optical fibers used, including 
the ones in the fiber array, are single mode fibers (SM800, 
core size of 5.6 um) designed for the wavelength range at 
around 780 nm. 


III. EXPERIMENTAL RESULTS 


The measured optical transmission spectra of the de¬ 
vice (outer radius of 10 um and ring width of 2.5 um) 


in air and water, which includes the insertion loss of the 
two input-output grating couplers, are shown in Fig. [3] 
(a). The peak transmission reaches ~ 0.3%, correspond¬ 
ing to an insertion loss of —13 dB per grating coupler. 
The coupling efficiency is limited by the separation be¬ 
tween the fiber probe and the grating couplers set by the 
thickness of the cover glass, which is around 100 um. In 
the measured spectra, groups of resonances having simi¬ 
lar quality factors and extinction ratios can be identified. 
They correspond to the TE whispering gallery modes of 
different radial orders. The resonance wavelengths are 
plotted against the coupling gap g in Fig. |3](b). Reso¬ 
nances of the same radial mode are plotted in the same 
color. Four modes can be identified in the spectra taken 
in air while three modes show up in case of water, which 
agrees with the simulation that the fourth TE mode in 
water is not well confined. 

Different radial modes show distinct behavior in qual¬ 
ity factor and extinction ratio. The loaded quality fac¬ 
tor Q and the normalized on-resonance transmission 
To = P 0 ut(A = Ao )/Pin are related to the intrinsic quality 
factor Qi and coupling quality factor Q c by 


1-1 = Q7 1 + Qc 1 

(1) 

_ (Q c - QA 2 

0 V Qc + Qi J 

(2) 


From the measured Q and To, Qi and Q c can be cal¬ 
culated, i.e. Q c ,Qi = Q± = ‘ZQ/Q- ± v^o)- From this 
it follows that <5- > Q+ but it cannot be determined 
from Eqs. (jT]) and (|2j) which one corresponds to Q c and 
which to Qi. This depends on whether the system is over¬ 
coupled ( Q c < Qi) or undercoupled ( Q c > Qi). One way 
to resolve this ambiguity is to observe how Q + and Q_ 
vary with the coupling gap g , since Q c is expected to 
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FIG. 3. (a) Optical transmission spectra taken in air and water, with coupling gaps of 100 nm and 310 nm respectively, (b) 

Resonance wavelengths of the devices plotted against coupling gap. Resonances belong to the same radial mode are plotted in 
same color, (c) Quality factors Q± of different radial modes plotted against coupling gap. 


be an increasing function of g while Qi is expected to 
be independent of g. Fig. [3] (c) plots the measured Q± 
against the coupling gap for different radial modes. For 
each radial mode, we attribute the Q± that has a weaker 
(stronger) dependence on g as the intrinsic (coupling) 
quality factor ( Q c ). Note that while Qi is expected to 
show no dependence on g , some of them do decay slowly 
as g decreases, which could be due to the increasing scat¬ 
tering loss in the proximity of the coupling waveguide. 
Also note that for the first mode in water, Q + and 
display a crossing feature: the system changes from the 
undercoupled regime at large g to being overcouplcd at 
small g. Among all the measured results, the highest 
loaded quality factor obtained is (1.53 ± 0.04) x 10 6 in 
air and (1.50 ± 0.03) x 10 6 in water. It is thus indeed 
possible to maintain the high Q that was obtained in air 
when operating the device in water. 

For an optomechanical device with such high optical 
Q , the optomechanical dynamical backaction can have a 
significant effect on the mechanical motion; it can either 
amplify or dampen the mechanical motion depending on 
the laser detuning l3fil . When the input laser power 
is high enough, the amplification effect can completely 
compensate the intrinsic mechanical damping and cause 
the resonator to self-oscillate, which has been shown to 
be useful in sensing application m ■ Because of the 
very high optical Q , the present device can be oper¬ 
ated in self-oscillation mode in air with low threshold 
power. Fig. [4] (a) shows the normalized transmission 


and the RF spectrum when the wavelength is scanned 
across the resonance with an input power of 45 uW. The 
transmission curve shows an asymmetric shape caused 
by the thermo-optical bistability [2E, where the heating 
due to the optical energy in the cavity moves the reso¬ 
nance to longer wavelengths. The transmission at even 
lower power (1 uW) is shown in the same figure for com¬ 
parison between the linear and nonlinear regimes. As 
the laser wavelength approaches from the shorter wave¬ 
length side of the resonance, the mechanical resonance 
peak becomes sharper and higher, indicating the effect 
of optomechanical amplification. Assuming the optical 
linewidth remains the same throughout the wavelength 
scan, the actual detuning can be inferred from the trans¬ 
mission value and the original linewidth. The mechani¬ 
cal frequency and the linewidth at different detuning is 
plotted in Fig. [4] (b). As expected, when the detuning 
approaches half the linewidth of the cavity (indicated 
by the dashed blue line), the mechanical linewidth de¬ 
creases until it reaches zero, where the device undergoes 
large amplitude self-sustained oscillation (region shaded 
in gray color in the figure) [39j . As the detuning is fur¬ 
ther reduced, the optomechanical effect diminishes and 
the mechanical linewidth gradually return to its original 
values until the thermo-optical bistable point is reached. 

When the device is immersed in water, the viscous 
damping is so great that we were not able to set the de¬ 
vice into self-sustained optomechanical oscillation. Nev¬ 
ertheless, with the very high optical Q, the thermome- 
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FIG. 4. Optomechanical oscillations in air. (a) Upper panel 
shows the normalized transmission across a resonance at dif¬ 
ferent laser power. Lower panel shows the contour plot of 
the measured RF spectrum versus wavelength when the laser 
power is 45 uW (corresponds to the black line in the upper 
panel). Cross-sectional plots of the spectra at different wave¬ 
lengths indicated by the dashed-lines are shown in the in¬ 
set. (b) Mechanical frequency and linewidth plotted against 
detuning from the optical resonance. The region where the 
device undergoes self-sustained oscillations is shaded in grey 
color. Blue dashed-line indicates the detuning that corre¬ 
sponds to the half-linewidth of the optical resonance. Black 
dotted-line shows the level of zero linewidth. 


chanical noise of the mechanical mode can be resolved 
even when the mechanical motion is highly dampened 
in water. Fig. [5] (a) shows the thermomechanical noise 
spectra of the mechanical resonance measured in both 
air and water environment. In both cases the thermal 
motion is visible as a peak on top of the detector impre¬ 
cision noise. In air, the peak is narrow and high; there 
the radial breathing mode has a resonance frequency of 
179.9 MHz and mechanical quality factor of Qm = 2160. 
In water, the frequency is slightly shifted down to 169.4 
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FIG. 5. (a) Thermomechanical noise spectrum of the radial 

breathing mode in air and water environment. Red solid line 
represents the fitting to the water data. Black and red dashed 
lines indicate the noise floors of the measurement in air and 
water respectively, (b) Resonance frequency (upper panel) 
and quality factor (lower panel) plotted against the width of 
the micro-wheel. r ou t is fixed as 10 um. The experimental 
(numerical) results are shown in black (red). 


MHz and Qm is reduced to 9. Note that the noise floors 
of the spectra taken in air and water are different, as in¬ 
dicated by the dashed lines in the figure, because a lower 
optical power was used in air to minimize the effect of 
the optomechanical backaction. For the noise spectrum 
in water, a noise floor of Shi? = 15 am/\/Hz is achieved. 

It is noteworthy that the resonance frequency drops 
by only 5.8% after immersing the device in water. This 
is very small compared to other demonstrated mechan¬ 
ical systems. For example, AFM cantilevers working at 
several hundreds of kHz can have resonance frequency 
that drop to 1/2 - 1/5 of the original (air) value when 
immersing in water m- Also silicon nitride nanostring 
resonators working at 100 MHz range show a more than 
25% drop in resonance frequency m- Another observa¬ 
tion is that, the quality factor of the present device is 
also higher than other demonstrated systems, which typ¬ 
ically have Qm of 5 or lower [9] fl5UT8i B0 I. The smaller 
frequency shift and higher Qm hints that the hydrody¬ 
namic mass loading and damping due to the surrounding 
water is relatively small for the micro-wheel resonator 
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structure. 

The dynamics of the micro-wheel resonator inside a 
liquid environment can be modeled with the equation of 
motion for a single harmonic oscillator after integrating 
over the mode shape: 

x(t ) + 2')±{t) + Fl 2 x(t) = Fftt) _|_ (t) , 

m m 

where f 2 , 7 and m are the resonant frequency, damp¬ 
ing rate and effective mass in vacuum, and Ff(t ) and 
Fth(t ) are the fluidic force and thermal fluctuation force 
acting on the resonator. In the limit of small oscilla¬ 
tion amplitude, the fluidic force can be assumed as a 
linear response function of the resonator motion, i.e., 
F f(t) = J dt'x(t - t’) x[t'). In the frequency domain, 
it can be written as Ff[u)\ = muj 2 F[uj\x[(jo\, where r[w] = 
—ix[w]/mw is a dimensionless function called the “hydro- 
dynamic function”. The definition adopted here differs 
from those in Refs. PH PS] by a geometrical factor. In 
general r[w] is a complex function. Its real part T^fw] 
and imaginary part T/fcu] are related to the mass load¬ 
ing and damping respectively due to the presence of the 
fluid. Using the fluctuation-dissipation theorem, it can 
be shown that the single-sided noise power spectral den¬ 
sity of the displacement due to thermomechanical fluctu¬ 
ation is given by pH |25] 

(AksT/mn^ur^u)} 

xx[1 (i-w 2 (i + r fl M))2 + (d; 2 r 7 H)2 5 u 

where Cj = uj/FL is the normalized frequency. Here the 
intrinsic damping rate 7 is assumed to be negligible com¬ 
pared with the hydrodynamic damping, i.e. 7 <C wT;[w], 

For a cantilever with rectangular geometry, the hydro- 
dynamic function can be solved numerically using the 
boundary integral method described in Refs. [2B10T1 PI2]. 
A formula introduced in Ref. [21] provides a good ap¬ 
proximation for the hydrodynamic function of a can¬ 
tilever with large width-to-thickness ratio. These studies 
are relevant to mechanical resonators with rectangular 
geometry such as AFM cantilevers. To study the hydro- 
dynamic effect on a micro-wheel resonator, here we use 
a boundary integral method for systems with rotational 
symmetry and solve the corresponding Green’s equation. 
The details of this model can be found in the supplemen¬ 
tary material. 

From the hydrodynamic function, the loaded resonance 
frequency f 2 W = H (1 + and quality factor 

Qw = (1 + r fl .[H. u; ])/r/[H u ,] in water can be calculated 
[2S] , Fig. [5] (b) plots the resonance frequency and qual¬ 
ity factor as function of ring width for devices with fixed 
outer radius r out = 10 um. The symbols are mean val¬ 
ues and the error bars are standard deviations for four 
devices per data point. As shown in the figure, the nu¬ 
merically calculated Qm is around 10 while the measured 
Qm is around 8.5 to 9, which agree fairly well with each 
other. The numerical mode also predicts the observed 
width dependence. 


The total loaded mass of the device including the 
water entrained to the resonator motion is given by 
m w = m(l + = 75 pg, which is roughly equal 

to the sum of the original (effective) mass of 65 pg and 
the mass of water enclosed in the Stokes boundary layer 
around the device (10 pg). The Stokes boundary layer is 
a measure of how far the fluid oscillations extend from the 
device and has a thickness 5 = 31 nm (See supplemen¬ 
tary materials). We would like to emphasize that the en¬ 
trained mass of water is only of a small fraction (15%) of 
the original resonator mass. This is orders of magnitude 
smaller than that of a typical cantilever structure, where 
the added water mass can be 10s of times higher than the 
resonator mass [19] 00]. Even for the very lightweight sil¬ 
icon nitride nanostring resonator, the added water mass 
is more than 1.3 times the original mass [151 . The small 
water entrainment of our device is the result of the small 
resonator dimensions, high resonance frequency and thus 
small Stokes boundary layer thickness, and in-plane mo¬ 
tion of the radial breathing mode, which is expected to 
have smaller mass loading effect in liquid compared to 
out-of-plane motion [12] , 

The values of both the resonator mass and entrained 
water mass are low, indicating that our device will have 
a large frequency shift A/ = ( f w /2m w )Am for a small 
added mass Am. This frequency shift is measured by 
driving the resonator to a coherent oscillation with am¬ 
plitude Xq and detecting the phase between the drive and 
the motion. The phase noise in such a measurement de¬ 
termines the mass sensitivity of the device 5m. The phase 
noise contains both the imprecision noise ( S nn ) and the 
resonator’s thermal motion ( S x x ). Both contributions 
can be estimated from the measured noise spectrum (cf. 
[ 5 ] (a)). The mass sensitivity of the resonator is given by 


(5) 


where the measurement bandwidth has been assumed 
to be much less than the mechanical linewidtli, i.e., 
Af l Fl w /Q w . The factor of 1/2 inside the square root 
accounts for the fact that only half of the noise power is 
distributed in the phase quadrature. From Fig. [5] (a), the 
sum of the thermomechanical noise and the measurement 
noise background gives a noise peak of 42 am/VHz on 
resonance. Assuming the device can be driven to an am¬ 
plitude of Xo = 100 pm (mechanical strain of 10 -5 ) before 
reaching the mechanical nonlinearity limit, the mass sen¬ 
sitivity is estimated to be 5m/(AFl/2Tr) = 2.5 ag/-\/Hz. 
This high mass sensitivity in water environment comes as 
a result of small mass loading and high quality factor of 
the device. So far, attogram sensing in fluid has only been 
achieved using cantilevers with embedded channels [13] . 
Our results suggest that optomechanical micro-wheel res¬ 
onator can be a promising candidate for in-situ attogram 
sensing in water environment. 


5m = 


m v 


Qw^O 


~ (^EX [f^TT;] T $7171) , , 

Z Z7T 
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IV. SUPPLEMENTARY MATERIALS 

A. Device fabrication 

This section presents the details of the fabrication pro¬ 
cess for the optomechanical micro-wheel resonator. A 
schematic illustration of the process flow is shown in 
Fig-! The device is patterned using the two-steps e- 
beam lithography and dry etching technique employed 
previously by our group [55] . The process starts with a 
silicon wafer with 2 um thermal oxide, 200 nm LPCVD 
stoichiometric silicon nitride, and a 100 nm thick PECVD 
silicon dioxide. First, e-beam lithography is performed to 
define the device pattern in an e-beam resist (ZEP520A). 
After developing, the e-beam resist is reflowed at 150 °C 
for 1.5 min to smoothen the sidewall. Reactive-ion etch¬ 
ing (RIE) with fluorine chemistry is then used to dry etch 
the top PECVD oxide and half (100 nm) of the silicon 
nitride layer. The e-beam resist is removed by O 2 plasma 
afterward and another ZEP520A resist is spin-coated on 
the chip for the second e-beam lithography. The second 
e-beam lithography defines the region where the bottom 
silicon oxide is to be exposed to HF solution for releas¬ 
ing. After the second lithography, another RIE is carried 
out to etch away the remaining half of the silicon nitride 
layer. During the second RIE, the top PECVD oxide 
functions as masking layer to protect the part of pho¬ 
tonic structure that is exposed to the plasma making the 
two patterns self aligned. Again, the e-beam resist is re¬ 
moved by O 2 plasma and the chip is further cleaned in 
piranha solution. 

PECVD SiO, #0 cover glass 

LPCVD Si,fJ 4 1 '' ^ 

SiO, 

I s ' | PDMS coating 

| "Two-step" lithography 

_ a _^ =aB=1 _ _ I Sand blasting 



| PDMS/Si0 2 bonding 

PECVD 



| Si Si0 2 Si 3 N 4 | PR Glass | PDMS tape 

FIG. 6. Fabrication process flow of the optomechanical 
micro-wheel resonator. The final result is highlighted by the 
orange box. 


Next, 3 um thick PECVD silicon dioxide is deposited 
on top of the chip. The chip is again cleaned in piranha 
solution and thoroughly dehydrated in an oven. It is 
then primed with HDMS and spin-coated with photore¬ 
sist (Shipley S1813). Photolithography is performed to 
pattern the microfluidic channels. The chip is post-baked 
at 120 °C for 10 min to further harden the resist. The 
chip is then immersed in a buffered HF solution. In this 
way, the etching of the microfluidic channels and the re¬ 
leasing of the micro-wheel is done in one step. Finally, 
the chip is transferred to isopropanol and flash dried on 
hotplate. At this stage, we have a chip of optomechanical 
devices located inside microfluidic channels (See Fig. [ 6 ] 
bottom right, and Fig. [l] (d) in the main text). 

The next step is to seal the microfluidic channels by 
bonding the chip to a cover glass slide. In order to mini¬ 
mize the distance between the fiber array and the on-chip 
grating couplers, a size #0 cover glass (80 - 100 um thick) 
is used, which is the thinnest among the standard glass 
slides. To bond the cover glass to the PECVD dioxide 
layer of the chip, a thin PDMS coating is used for adhe¬ 
sion. PDMS is prepared by mixing the PDMS elastomer 
and the curing agent in 10:1 ratio. The mixture is then 
diluted in hexane with 1:1 vol. ratio and is thoroughly 
mixed using a vortex mixer. The diluted PDMS solution 
is spin-coated on a pre-cleaned cover glass at 3000 rpm 
for 30 s and is cured overnight in an oven at 100 °C. 
This results in a 300 nm thick film. After the PDMS is 
completely cured, an adhesive tape punched with holes 
(diameter of 2 mm) was taped on the glass side to act as 
the mask for the subsequent sand blasting process. Sand 
blasting with AI 2 O 3 powder is used to drill holes in the 
cover glass. The dust generated during the process can 
be removed using adhesive tape. To proceed to the bond¬ 
ing, the glass piece and the chip are put in O 2 plasma 
with an RF power of 50 W for 30 s. The two pieces are 
then aligned and pressed against each other. The bonded 
chip is put in the oven for overnight baking. For the final 
packaging, tubing adapters are glued to the cover glass 
at the hole area using epoxy. Tubings are attached to the 
adapter for liquid transfer. 


B. Wavelength calibration 

To quantify ultra-high Q optical resonances, precise 
calibration of the laser wavelength is very important. The 
tunable diode laser used in the measurement (New Focus 
TLB-6712) has a piezo tuning mechanism to fine tune the 
cavity length and thus the laser wavelength. To calibrate 
the exact wavelength shift induced by the voltage applied 
to the piezo, we use a wavelength meter (HP86120) to 
measure the laser wavelength. Fig. [7] shows a calibra¬ 
tion curve of the laser wavelength as the piezo voltage 
is varied. Fluctuation in the wavelength readings is due 
to instrument noise, which depends on the measurement 
bandwidth. The result shows that the laser wavelength 
can be tuned with a high degree of linearity and precision. 
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Piezo tuning (%) 

FIG. 7. Laser wavelength plotted as a function of the piezo 
voltage (expressed as a percentage of the full range). 


C. Numerical model of hydrodynamical loading 

This section describes the theoretical model for calcu¬ 
lating the hydrodynamic function of systems with rota¬ 
tional symmetry. The hydrodynamic effects on an oscil¬ 
lating solid structure can be understood using the follow¬ 
ing physical picture. As the boundary of the structure 
oscillates with a certain amplitude A , the surrounding 
fluid is displaced and the motion propagates through the 
medium with dynamics governed by the Navier-Stokes 
equations. Eventually the whole system reaches a steady 
state. The fluid force acting on the solid structure can be 
calculated from the pressure field and the viscous drag¬ 
ging of the fluid. When the oscillation amplitude is small, 
the fluidic force is expected to be linear in x(t) and, in 
general, to have a relative phase shift compared to the 
mechanical motion. The phase shift means that part of 
the force is in-phase with the velocity which causes damp¬ 
ing, and part of the force is in-phase with the displace¬ 
ment which causes mass loading. The above scenario 
can be formulated as a partial differential equation prob¬ 
lem of solving Navier-Stokes equation. For the boundary 
condition at the fluid-structure interface, the no-slipping 
condition is assumed, which means that the velocity of 
the fluid at the surface follows the velocity of the me¬ 
chanical motion of the solid structure exactly. We also 
assume the fluid to be incompressibile, i.e., V • v = 0, 
where v is the velocity field of the fluid. 

Reynold’s number is an important scaling parameter in 
fluid dynamics. In general, one can define two Reynold’s 
numbers based on the ratios between the inertial force 
pjyiV = p-§iv + p{v ■ V)fT and the viscous force /iV 2 u. 
For oscillatory motion, the time-dependent term of the 
inertial force gives rise to the frequency-based Reynold’s 
number Re/ = puiL 2 / p and the nonlinear term of the 
inertial force gives rise to the velocity-based Reynold’s 
number Re u = pujAL/p [14], where p and p are the den¬ 
sity and dynamic viscosity of the fluid, L is the char¬ 
acteristic length scale of the system, and A is the os¬ 
cillation amplitude. One can immediately see that for 



FIG. 8. (a) Cross-sectional view of the micro-wheel structure 
in cylindrical coordinates, (b) Real and imaginary parts of the 
calculated hydrodynamic function plotted against frequency. 

small-amplitude oscillations, i.e., A/L <C 1, the nonlin¬ 
ear term becomes negligible and therefore the linearized 
Navier-Stokes equation can be used. As an example, 
the magnitude of the thermo-meclranical fluctuations is 
A ss (x 2 ) ~ 0.1 pm while the characteristic length is 
L ~ 1 pm (see below). 

The linearized Navier-Stokes equation in the frequency 
domain can be written as 

- ^v = --VP + V 2 u. (6) 

o z n 

Here, 5 = \Jp/pu is the Stokes boundary layer thickness, 
which is related to the frequency-based Reynold’s number 
by Re/ = ( L/S ) 2 . Following the convention of Refs. [14, 
1213; L = {rout — ri n )/2 is chosen as the half-width of 
the micro-wheel. For the device with width of r out — 
ri n = 2.5 pm and frequency of 170 MHz, it gives S = 
31 nm and Re/ = 1700, which is very high compared 
with other demonstrated systems. Note that the velocity- 
based Reynold’s number Re M is still much less than one. 
Therefore, here we have Re u <C 1 <C Re/. 

Consider the system shown in Fig. [8] (a). All mechan¬ 
ical motion is assumed to be rotationally symmetric, i.e., 
v(r) = v r {r,z)f + v z (r,z)z. The velocity of the micro¬ 
wheel structure is given by e~ lult v s f while the velocity 
field of the surrounding fluid is given by e~ lut v{r , z ). 
The no-slipping condition implies that v(r. z) = v s f 
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for (r, z) £ C at the boundary. For systems that ex¬ 
hibit rotational symmetry, it is convenient to use the 
stream function ip to describe the fluidic motion, where 
V x ip = v. Another useful quantity is the vorticity de¬ 
fined as w = V x v. It can be shown that both ip and 
w have only azimuthal component, i.e., ip = ip(r, z)9 and 
w = w(r,z)8. By rewriting the Navier-Stokes equation 
in terms of ip and w and applying Green’s identity to 
convert the differential equation into integral equation, 
it can be shown that ip satisfies the following equation 


iP(r\z') 


= <b \dl\ 


c 


,P. 


ip ( d n G +^f-n) - Hd n ip 


$ 


—i6 —dt$ + id w ( d n $ -\—f • n 


(7) 


where d n and d t represent the differentiation along the 
normal and transverse direction. The Green’s functions 
G, H, $ are given by 


G{r',z'\r,z) = - 
7r , 

H(r', z'\r, z) = - 
7r t 

z'\r, z) = G( 


nOO 

/ dk cos k(z — z')Ii(kr < )K 1 (kr > ) 

Jo 

dkcosk{z — z')Ii(kr < )Ki(kr > ) 

( 8 ) 


r',z'\r, z) - H(r', z'\r, z) 


where I\(x) and Ki[x) are the modified Bessel function of 
the first and second kind with index 1, k = ^k 2 — i/8 2 , 
r< (r>) represents r or r' whichever is smaller (large), 
and h is the unit vector point outward from the surface 
(see Fig. [8] (a)). For computation of the above Green’s 
functions with large arguments, the following asymptotic 


form of Ii(x)Ki(x') becomes useful 


Ii{x)K\(x') 


1 A A 

8 a: 8x' 


15 


15 


128a; 2 128a/ 2 


9 

64a; x' 



(9) 


since it ensures that the product exponentially decays 
with x and x' as long as x' > x and thereby avoids the 
overflow problem when computing I\ with large argu¬ 
ments. 

Eq. (7| means that the stream function ip anywhere 
inside the fluid is fully determined by the values of ip, 
d n tp, P, and w at the boundary C. Among these four 
functions, ip and d n ip are known since they can be cal¬ 
culated from the velocity of the fluid at the boundary, 
which follows that of the mechanical resonator because 
of the no-slipping boundary condition. When the term on 
the left hand side evaluated at the boundary surface, i.e., 
(r',z')GC, Eq. 0 represents a self-consistent equation 
with two unknown functions P(r, z) and w(r, z). Another 
independent equation can be obtained by taking gradi¬ 
ent d n / on both sides of the equation. The two integral 
equations with two unknown functions can be solved, for 
example, using the discretiziation method presented in 
Refs. ?Ti~. Ti l . 

After solving P(r,z) and w(r,z) at the boundary, the 
fluidic force acting on the structure surface can be calcu¬ 
lated from the equation m, 


Ff = j> (—Pn + nw x h)dA . (10) 

The total fluidic force is a sum of the pressure force 
acting perpendicular to the surface and the viscous 
force acting in parallel to the surface. From the flu¬ 
idic force the hydrodynamic function can be obtained, 
i.e. Ff[co] = mw 2 r[wji[w]. Fig. [5] (b) plots the calcu¬ 
lated hydrodynamic function in the same frequency range 
as the measured data in Fig. [5] (a). Water density of 
p = 1000 kg/m 3 and viscosity of p = 1 mPa • s are used. 
The device dimension are r ou t = 10 pm and Vi n = 7.5 pm 
in this calculation. 













